Model Independent Early Expansion History and Dark Energy 



Johan Samsing 1 ' 2 , Eric V. Linder 1,3 , Tristan L. Smith 1 
1 Berkeley Center for Cosmological Physics & Berkeley Lab, 
University of California, Berkeley, CA 94720, USA 
2 Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, 
Juliane Maries Vej 30, 2100 Copenhagen, Denmark 
^Institute for the Early Universe WCU, Ewha Womans University, Seoul, Korea 

(Dated: August 27, 2012) 

We examine model independent constraints on the high redshift and prerecombination expansion 
history from cosmic microwave background observations, using a combination of principal component 
analysis and other techniques. This can be translated to model independent limits on early dark 
energy and the number of relativistic species N e s- Models such as scaling (Doran- Robbers), dark 
radiation (AJV e fi), and barotropic aether fall into distinct regions of eigenspace and can be easily 
distinguished from each other. Incoming CMB data will map the expansion history from z = 0-10 5 , 
achieving subpercent precision around recombination, and enable determination of the amount of 
early dark energy and valuable guidance to its nature. 



I. INTRODUCTION 

The expansion history of the universe is a fundamen- 
tal property of cosmology, reflecting the energy density 
constituents and their evolution. Yet remarkably little is 
known in detail about it, other than in a coarse grained 
average. For redshifts between 3000 and 10 9 , the uni- 
verse was mostly radiation dominated, for redshifts be- 
tween 3000 and ~ 1 it was mostly matter dominated, but 
excursions are possible - in the effective number of rela- 
tivistic species N e & say, or even temporary breakdown of 
such domination - and the level of subdominant compo- 
nents is not well constrained. Only around the epoch of 
primordial nucleosynthesis and of recombination is the 
expansion rate (Hubble parameter) better constrained, 
but even there at the ~ 5% level averaged over the epoch 
[1.2]. 

Given the importance of the expansion history, and the 
improvement in cosmic microwave background (CMB) 
data, we investigate what constraints can be placed on 
it in a model independent way, i.e. other than fitting for 
a deviation of a particular functional form such as extra 
7V e ff or a specific dark energy model. This would fill in a 
vast range of cosmic expansion where almost no precise 
constraints have been placed. That is, an error band for 
the Hubble parameter H(z) at z > 1000 should be a 
staple of cosmology textbooks, and yet does not exist. 

The early expansion history has an important bear- 
ing on understanding the nature of dark energy as well, 
the question of 'persistence of dark energy. For a cosmo- 
logical constant A, the dark energy density contributed 
at recombination is J7a ~ 10 -9 , while the current up- 
per limit from data is above 10 -2 . This gives substantial 
unexplored territory. Moreover, the current constraints 
use a specific functional form for the dark energy evo- 
lution (usually the Doran- Robbers form [3]), but other 
models could lead to significantly different limits [4]. 
Thus, model independent limits on early dark energy are 
needed. Physics origins for early dark energy can be quite 
diverse, e.g. from dilaton models (as in some string the- 



ories) to k-essence (noncanonical kinetic field theories) 
to dark radiation (as in some higher dimension theories) 
[5]. Establishing whether CMB observations could dis- 
tinguish these classes is another important question. 

Improvement of CMB data recently by higher resolu- 
tion observations extending the temperature power spec- 
trum to multipoles I ~ 3000 by the Atacama Cosmology 
Telescope (ACT [6]) and South Pole Telescope (SPT [7]) 
gives valuable leverage since higher multipoles are sensi- 
tive to modes crossing the cosmological horizon at earlier 
times. This advance was used in [8] to rule out in a model 
independent manner the presence of any epoch of cos- 
mic acceleration between z ss 2 and 10 5 (supplementing 
the limits from growth of structure post-recombination 
in [9]). Upcoming Planck and ground based polarization 
experiment data will also map out the polarization power 
spectra, giving additional constraints. 

To carry out a model independent analysis of the early 
expansion history, we use a combination of redshift bin- 
ning and principal component analysis. In Sec. II we 
lay out the methodology for describing arbitrary H(z). 
Analyzing the results in Sec. Ill, wc identify the redshifts 
ranges where the CMB observations are most sensitive to 
expansion variations. We project three classes of mod- 
els representing different physical origins onto the eigen- 
modes to explore the discriminating power of the data 
in Sec. IV. In Sec. V we discuss the results and future 
prospects. 

II. EXPANSION HISTORY 

The expansion Friedmann equation directly relates the 
expansion rate of the universe, or Hubble parameter, to 
the energy density constituents, 

where we neglect curvature (from a theoretical prior for 
flatness and because we mostly treat high redshift where 



2 



it would be negligible). At high redshift the canonical 
expectation is that the universe is matter or radiation 
dominated, so we write 
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Deviations 5(a) to the fiducial expansion rate can also 
be interpreted as an effective dark energy density differ- 
ing from that of the cosmological constant, with 
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We can write the dark energy density evolution as 
Pde(a) = pdefi f(a) = pa f(a) 
f(a) = l+(l+^\s, 
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where pb g is the background energy density excepting 
dark energy, i.e. usually the dominant component, matter 
or radiation. More simply, the fraction of critical density 
contributed by the effective dark energy is 
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We can readily see that at high redshift we obtain a frac- 
tional early dark energy density contribution of approxi- 
mately 5(a), for f^A(a) <C 5 <C 1. During epochs when 5 
is constant, this is a constant fractional contribution. 

Our goal is to analyze constraints on the variations 
5(a) from the canonical model with 5 — 0. We begin by 
writing 5(a) as a linear combination in an orthogonal bin 
basis, 



5(a)=Y J PMa) , 
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where bi(a) is a tophat of amplitude 1 over a given range 
of scale factor a, and otherwise. That is, the Hubble 
parameter deviations 5(a) are given as a linear combina- 
tion of piecewise constant values. We can then constrain 
H(a) in bins of a, a model independent description. The 
bin basis is also the standard first step in principal com- 
ponent analysis (see, e.g., [10]), as we will pursue in the 
next section. We choose N bins per decade of scale factor 
over the range of log a = [—5, 0], beginning with N = 20. 

Since we are interested in the expansion history we deal 
directly with the Hubble parameter (or effective dark en- 
ergy density). Treating bins of the dark energy equation 
of state, or pressure to density, ratio w(z) would have 
some drawbacks here. Most severe is that to obtain H(z) 
one must integrate w(z') over all redshifts from zero to 
z. This makes it difficult to explore the early expan- 
sion history in a model independent manner. Moreover, 
the instantaneous w(z) is not fully informative: during 
matter domination, for example, any level of dark energy 



density from Vtd e = 10~ 9 to 10 -2 or whatever that scales 
as the matter has w — 0. Thus we aim to derive con- 
straints directly on variations in H(z), and consider the 
interpretation of these as a further step. 

The expansion history directly feeds into the CMB 
power spectra, through changing the distance scales, e.g. 
of the sound horizon or damping scale, and the rela- 
tion of multipole i (or angular scale 9) to wavenumber 
k — £/rj(z), where r\ is the conformal distance. It also 
affects the decoupling of photons from baryons and the 
growth of perturbations in both. 

The treatment of perturbations requires some atten- 
tion. The description of the cosmic expansion gives the 
evolution of the homogeneous background, but consis- 
tency of the field equations requires consideration of per- 
turbations in all components of energy density. Unless 
the deviation in H(z) is interpreted purely in terms of 
a cosmological constant (which indeed is purely homoge- 
neous), spatial perturbations have to be accounted for, at 
least formally and generally in actual practice. The per- 
turbation evolution equations for the additional energy 
density involve the quantities w(z), w'(z) — dw/dlna, 
the initial conditions on the density perturbation, and 
the sound speed of the effective fluid c s (z). (One could 
also add a viscous sound speed or anisotropic stresses, 
see [11].) 

The first three of these are fairly straightforward. For 
any deviation 5(a) one can define an effective equation of 
state 
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where the last line holds when p\ gives negligible contri- 
bution to the effective dark energy density. One can take 
a further derivative to obtain w' . Initial conditions are 
usually taken as adiabatic and stresses are taken to van- 
ish. However, one does have to specify the sound speed. 
If one interprets the extra energy density as arising from 
quintessence, i.e. a minimally coupled, canonical scalar 
field, then c s — 1. In general, the necessary inclusion of 
perturbations in whatever is the origin of the deviations 
in the expansion history prevents a purely model inde- 
pendent treatment - one has to make some assumptions 
about the physics. Here we fix c s to that for the partic- 
ular cases we consider, but in future work we will fit for 
it. 

We modified CAMB [12] to allow a general H(z), with 
the w(z) that goes along with this. We then solve the cou- 
pled background evolution, and photon, matter, and ef- 
fective dark energy perturbations equations to obtain the 
CMB power spectra. For evaluating binned H(z) models, 
using the orthogonal bin basis introduced in Eq. (9), we 
slightly smooth the bin edges, using a Gaussian smooth- 
ing of width 0.075 times the bin width, to prevent in- 
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finite derivatives. We extensively test convergence and 
numerical stability of the results (also see [8] where this 
procedure was found to be robust). 

Figure 1 shows the bins in log scale factor (20 bins 
per decade) and the conversion to multipole space (over- 
laid with the CMB temperature power spectrum) by I ps 
770/77(0), which approximately relates a given wavenum- 
ber k s» £/tjq to the time it entered the horizon. Note that 
a uniform binning in log a, which is the expected charac- 
teristic scale for physical variations in the expansion, is 
not uniform in multipole space. 

To place constraints on allowed deviations 5(a) we 
carry out a Fisher matrix calculation. The Fisher ma- 
trix elements are given by 

where X, Y is any combination of the CMB tempera- 
ture power spectrum (T), E-mode polarization power 
spectrum (E), and temperature-polarization cross power 
spectrum (TE). The covariance matrix COV is given 
by the measurement uncertainties of the CMB observa- 
tions; we adopt the characteristics of the Planck satel- 
lite experiment [13]. The parameter set {p^ includes the 
usual cosmological parameters - the physical baryon den- 
sity ftbh 2 , physical cold dark matter density £l c h 2 , total 
present matter density tt m (the present Hubble constant 
h is a derived quantity), primordial scalar perturbation 
power law index n s , optical depth r, and present ampli- 
tude of mass fluctuations as - and the iVbin expansion 
variation parameters 6(a.i). The uncertainties on each 
5(o x -) are given by the (square root of the) respective di- 
agonal element of the inverse of the Fisher matrix. 

Figure 2 shows the sensitivity y^M? °f the weighted 
combination of CMB power spectra (T, E, TE) to the 
expansion deviations in each redshift bin for each mul- 
tipole. Sensitivity peaks around the acoustic peaks and 
is reduced at low multipoles due to cosmic variance and 
at high multipoles due to the finite resolution from the 
instrument beam size. Since the polarization power spec- 
trum is out of phase with the temperature power spec- 
trum, dips in the temperature sensitivity are filled in by 
polarization information. 

Figure 3 shows the actual Fisher and covariance sub- 
matrices corresponding to the expansion bin parameters 
(marginalized over other parameters in the case of the 
covariance matrix). First, we notice the maximum of the 
information content is near decoupling (log a —3), as 
expected. Earlier times, log a < —4, map to multipoles 
on the damping tail and so have less leverage, while re- 
cent times, log a > —3, are on the Sachs- Wolfe plateau 
and again have limited information. The Fisher matrix 
is not diagonal because expansion deviations affect all 
later times, e.g. perturbation evolution once the wave- 
mode is within horizon and integral quantities such as 
the sound horizon. This will be one of the motivating fac- 
tors for carrying out principal component analysis (PCA) 
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FIG. 1. [Top] The bin basis for expansion variations 8(ai) is 
plotted vs log a. [Bottom] The scale factors Oj of the center of 
each bin in the top panel (with matching colors) are approx- 
imately mapped into multipoles by I ~ rio/rj(ai), with the 
CMB temperature power spectrum overplotted. (Amplitudes 
of the colored lines are arbitrary.) 



in Sec. III. The covariance matrix (inverse of the Fisher 
matrix), however, has a more diagonal structure, and so 
bins can be a useful parameter set if carefully chosen (see 
Sec. V). 

Most importantly, Figure 4 shows the constraints on 
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FIG. 2. The fractional Fisher sensitivities, »/ Faj, to the expansion variations S(ai) in each bin, color coded as in Fig. 1, are 
plotted vs multipole. The top panels use log scale in multipole, the bottom panels a linear scale to highlight different regions. 
The left panels are for the temperature power spectrum only, while the right panels use the variance weighted sum of the T, 
E, and TE power spectra entering the Fisher information matrix. Polarization information fills in the sensitivity gaps due to 
the acoustic troughs. 



the expansion history 

a{H) 



a{8) 
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i.e. the fractional uncertainty on H (a) due to deviations 
<5, marginalized over the other cosmological parameters. 



This is the "textbook" plot, showing the state of our 
knowledge of the early expansion history of the universe 
when given CMB data of Planck sensitivity. The con- 
straints depend on what bandwidth we wish to constrain 
the expansion history: the top curve shows 10 bins per 
decade, the bottom curve 2 bins per decade. One can 
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FIG. 3. The Fisher information submatrix (left) corresponding to the expansion variation parameters 5(a,i) is plotted with 
redder shades representing larger absolute values (more information). The color bar gives the log of the absolute value of 
elements. The covariance matrix (middle), marginalized over other cosmological parameters, follows the same color scheme, so 
the best determined parameters (smallest errors) are bluest. Diagonal elements of the Fisher matrix (right), as a fraction of 
the largest diagonal element, quantify for which redshifts the CMB data is most sensitive to the expansion history. The bump 
at log a w — 1.1 reflects reionization. 



trade off sensitivity to fine features vs overall level of con- 
straints. With 10 bins per decade one can achieve percent 
level constraints on H 2 (a) near decoupling, while with 2 
bins per decade one obtains subpercent constraints over 
more than two decades in scale factor. The relation be- 
tween 10 bins and 2 bins is not simply a %/5 scaling due to 
correlations between bins (the offdiagonal elements of the 
covariance matrix), and the lowest redshift bin is partic- 
ularly affected by covariance with the other cosmological 
parameters. 



III. PRINCIPAL COMPONENTS OF 
EXPANSION 

Expansion deviations at some redshifts may have sub- 
stantially the same consequence for the observables as 
deviations at some other redshift, or deviations may be 
correlated in such a way that only the difference between 
them is important. This leads to the idea of compressing 
the 100 bins between log a = [—5,0], or at least the in- 
formation contained in them, into fewer variables. One 
might for example speculate that the major effect of ex- 
pansion deviations for a > 10~ 3 comes from shifting the 
distance to CMB last scattering, and whether the varia- 
tion occurs at a = 0.01 or a — 0.1 is less crucial. 

Principal component analysis (PCA) can provide an 
efficient way to compress the influence on the observables. 
For some uses probing dark energy and the CMB, see for 
example [10, 14-18] (and dark matter and the CMB in 
[19]). By diagonalizing the Fisher matrix we can find its 
eigenvectors that best summarize the sensitivity of the 
observations to the expansion deviations. We can then 
transform the bin basis to an orthogonal eigenmode basis 
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FIG. 4. The fractional precision with which the expansion 
history can be determined by projected Planck CMB data is 
plotted vs scale factor, for two different bandwidths. The top 
(bottom) curve is for 10 (2) bins per decade in scale factor. 
Subpercent precision can be achieved around decoupling but 
large swaths of the cosmic history will remain unknown. 

of principal components (PCs), writing 

6(a) = J~]mi ej(a) , (15) 

where m, is the amplitude of mode i and e^(a) is the 
eigenvector. Since the modes are orthogonal, the errors 
(Tj = er(mj) on the amplitudes are uncorrelated. Using 
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the entire set of bins or the entire set of PCs is equivalent, 
but using only a few PCs with the highest eigenvalues 
(smallest cr,) in general allows one to approximate the 
full set more accurately than the same number of bins. 
That is, the information can be compressed. 

Figure 5 illustrates some of the PCs for the CMB 
Fisher matrix, ordered from highest to lowest eigenval- 
ues (best to worst determined). Note that as expected 
most of the activity in the first PCs is prerecombination, 
associated with the acoustic peaks. In modes 7-10, there 
is some low redshift action coming primarily from the in- 
tegrated distance to last scattering and the reionization 
epoch. Higher PC modes tend to be more oscillatory (es- 
sentially high derivatives of the expansion behavior) and 
localized. 

By taking the cumulative sum of the eigenvalues, we 
find that the first 7 PCs contain 99% of the variance. 
That is, Ei °"i~ 2 / Ei°° o"i" 2 = 0-99. This means that 
the great majority of expansion behaviors, as far as their 
observational detectability is concerned, can be described 
with just 7 parameters mi, . . . m-j. 

It is convenient to normalize the PCs such that 
En e i(°n) e j K) = where dij here is the Kronecker 
delta, and then the mode coefficients m,j give the ampli- 
tudes for a given model of expansion history deviation, 



= 22S(a n )ei(a n ) , 



(16) 



where a n denotes the bin centers. For many narrow bins 
the sum can be converted to an integral. The amplitudes 
TOj only have meaning when discussing a specific model; 
note the canonical model ACDM has all m ; = 0. Since 
the mode amplitude rrii has no a priori magnitude, we 
do not know whether its uncertainty Ui — 0.01, say, is 
a small or large number. We can compare <7j's to each 
other (but if the to,; are this is irrelevant), or for a 
specific model we can form a signal to noise ratio for a 
mode, 
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If (S/N)i -C 1, then it does not matter if <jj looks small, 
the mode cannot be well measured. Thus caution must be 
used in interpreting PCA (for more details see [17, 20]). 

Although PCA certainly compresses the information 
content into fewer parameters, the number of PCs re- 
quired to describe arbitrary expansion histories is still 
large for the purposes of, say, Monte Carlo simulations. 
We can therefore use PCA instead as a guide in two ways: 
to examine the ability to discriminate among different 
classes of models for the expansion deviation, and to in- 
dicate which regions in a show the most sensitivity to 
deviations and hence which of the original bins are most 
useful. These are treated respectively in Sec. IV and 
Sec. V. 
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FIG. 5. Principal components of the CMB observational sen- 
sitivity to expansion history are plotted vs log a, for the first 
10 and last 2 modes. The bottom plot shows the at for all 
modes. 
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IV. COMPARING EARLY DARK ENERGY 
MODELS 

While we wish to concentrate, as much as possible, on 
a model independent approach to constraining the ex- 
pansion history, it is useful to make contact with various 
classes of models to make sure that important behaviors 
are captured. In broad strokes, one can consider cases 
where the expansion deviations decline at times earlier 
than recombination, increase at earlier times, or remain 
constant. This can be translated into early dark energy 
models that contribute a lower fraction of the total en- 
ergy density in radiation vs matter domination, a higher 
fraction, or a constant fraction. This tilt of the expan- 
sion rate can be an important discriminant, somewhat 
analogous to the tilt of the primordial power spectrum 
for inflation. 

Early dark energy models have been proposed with 
each of these behaviors. Examples of the three classes 
respectively are 1) Barotropic dark energy with sound 
speed c 2 = [21], sometimes called aether models, 2) 
Barotropic dark energy with sound speed c 2 — 1/3 [21], 
sometimes called dark radiation, and 3) Scaling dark en- 
ergy such as the commonly used Doran-Robbers model 
with c 2 = 1 [3]. 

Barotropic models have w(a), pd e (a), and hence 5(a), 
all determined by c 2 s (a). The dynamics is given by [21] 



u/ = -3(1 + w)(c 2 - w) 

with solution for constant c s of 

, , c^a- 3 ( 1+c '' - 1 
W(fl) = Ba -3(i-H*) + i 



(18) 



(19) 



where B = (1 + w ) / (c 2 — w ) . For the c? s = aether 
model, 



Pae(a) = Poo + tt e p m . a 3 

x , s Q m (a" 3 -1) 
d ae (a) = il e 



(20) 
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Poo/Pcritfi = 1 - ^m(l + ^e) = (-™o)(l - ^m) (22) 
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(23) 



where during matter domination the dark energy con- 
tributed a constant fractional density Q e , but this de- 
clines at earlier times as radiation becomes important. 
For the c 2 — 1/3 dark radiation model, 

Pdr(a) — Poo + QePrfia^ 4 

a r ( a - 4 - 1) 

fi r (a~ 4 - 1) + n m (a- 3 - 1) + 1' 

Poo/ Pcritfl ~ 1 — f^m (26) 

wo « -1 , (27) 

where during radiation domination the dark energy con- 
tributed a constant fractional density Q e , but this de- 
clines at later times as matter becomes important. 



Sdr(a) = fi e 
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The most commonly used early dark energy is the 
Doran-Robbers form [3], 
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where during both matter and radiation domination the 
dark energy contributed a constant fractional density Q e . 
The sound speed is conventionally taken to be c 2 = 1. 

Thus the three models we consider have expansion his- 
tory deviations with complementary behaviors: rising, 
falling, and constant. This range should give a good 
indication of how PCA can characterize the expansion, 
while also having physical motivations. Of the physics 
origins mentioned in the Introduction, some string the- 
ories give Doran-Robbers behavior, some noncanonical 
kinetic fields give aether behavior, and some higher di- 
mension theories give dark radiation. 

Once we calculate the m, for a model through Eq. (16), 
we have a better indication of the importance of a PC 
mode through the signal to noise criterion of Eq. (17) - 
recall that the <7j alone say little about whether the mode 
is relevant. To choose the number of modes to keep in 
describing a model, we can simply impose a S/N cutoff, 
or ask that the cumulative S/N of the modes kept be 
above some fraction, e.g. 95% of the total S/N from all 
modes. 

Another indicator is the statistical risk, or mean 
squared error. This takes into account that more modes 
decrease the bias from the true model, but increase the 
accumulated variance in 5. The risk i?j\r is 



Rn — yb N 
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One could choose to keep that number N PCs where Rn 
is minimized. 

Figure 6 shows for each of the three models the approx- 
imation to 6(a) as more PCs are added, the values mi, 
and the cumulative S/N and risk as a function of number 
of PCs used. In general we find the risk requires more 
PCs than the S/N criterion; this makes sense since S/N 
concentrates on those PCs fitting the observable (CMB 
power spectra) while risk attempts to fit the unobserv- 
able expansion deviation. Another drawback to risk is 
that it does not scale with the amplitude of the modes, 
i.e. while S/N increases linearly with J7 e , the bias term 
in the risk scales but the variance does not, so the risk is 
weighted unevenly depending on deviation amplitude. 

For the dark radiation, aether, and Doran-Robbers 
models, respectively, we should keep the first 6, 7, and 5 
PCs according to S/N, and 10, 14, and 15 PCs according 
to risk. However, we note that this is if we keep the PCs 
in order according to <x;. If we choose the highest S/N 



modes individually, we only require 4, 3, and 3 modes to 
attain 95% of the full S/N (but this requires knowing the 
correct model, or assuming a given set of models). 

Let us examine these models in more detail. Figure 7 
shows how these models are well separated in eigenmode 
coefficient space. Considerable discriminating power oc- 
curs using modes 1 and 7, for example, with the sepa- 
rations between the three models many times the uncer- 
tainties cri. From the shape of the modes in Fig. 5, we can 
see that mode 1 is roughly measuring the amplitude of 
the expansion deviation at recombination, and whether it 
is increasing or decreasing (thus distinguishing all three 
models), and mode 7 is sensitive to more recent devia- 
tions such as occurring in Doran-Robber and aether, but 
not dark radiation, cases. Thus these two modes together 
are adept at distinguishing between the rising, falling, 
and constant deviation classes of expansion history, and 
early vs late deviations. 

We emphasize that the PCs are really telling us about 
fits to the observable power spectra and not reconstruc- 
tion of the expansion deviation in a fine grained sense. 
While Fig. 6 shows that ~ 50 PCs are needed to model 
5(a) well for these example, Fig. 8 demonstrates that 
only ~ 10 PCs are needed to give accurate agreement in 
Cg. Figure 8 shows the residual of the CMB tempera- 
ture power spectrum from the sum of the first N PCs 
(in S/N ordering, and also shown in er, ordering for the 
aether model) relative to the true model for the three 
models, as well as the x 2 summed over multipoles, and 
the reconstructed expansion history. Although the recon- 
structed expansion history may only agree over certain 
ranges of o, this can still give excellent agreement for the 
observables as the power spectra are not equally sensitive 
to all scale factors. 

The numbers of PCs required to ensure an accurate 
estimation of the observables, say x 2 < 10 (summed over 
3000 multipoles), is generally greater than 10, rendering 
cumbersome a straight application of mode amplitudes 
as parameters in a Monte Carlo simulation. (And recall 
that for a model independent analysis we do not know 
a priori the S/N ordering, so one would need to include 
many more PCs using ai ordering.) Moreover, PCs per 
se are not always easy to interpret in terms of physical 
effects. Therefore it is both clearer and more efficient 
to use the PCA instead to guide an effective, low order 
binning basis. 



V. CONCLUSIONS 

Our knowledge of the expansion history of our uni- 
verse, even at the level of degree of matter domination 
or radiation domination at early epochs, is remarkably 
imprecise. Cosmic microwave background radiation mea- 
surements from ACT, Planck, and SPT (and later ACT- 
pol and SPTpol) will shed light on the times around re- 
combination and reionization. We quantify the model in- 
dependent state of our knowledge through a combination 



of redshift bin and principal component analysis, finding 
that subpercent level constraints will be placed by Planck 
over log a = [—2.5, —5] for a bandwidth of A log a = 0.5. 

CMB data will address one of the key aspects of dark 
energy - its persistence, a characteristic of many high en- 
ergy physics origins - and we find that several different 
classes of early dark energy are well separated in PCA 
space. The limits can also be interpreted in terms of 
the number of effective relativistic species, iV c ff, such as 
an extra neutrino type, with current data mildly prefer- 
ring further contributions. A thermal relativistic neu- 
trino species adds 23% to the photon energy density, so 
6 = 0.13 AiVeg, giving tight limits on extra relativistic 
degrees of freedom from the forthcoming data. 

We explore three specific models, representing differ- 
ent classes for early dark energy, possibly corresponding 
to different physical origins. The commonly used Doran- 
Robbers form has a dark energy fraction f2 e that is con- 
stant through the recombination epoch. We also investi- 
gate a dark radiation model with 17^ (a) rising to the past 
and a barotropic aether model with f2<j e ( a ) falling to the 
past, and find that the dominant PC mode is well able 
to distinguish between these behaviors. Since the am- 
plitude of that mode is greatest for the Doran-Robbers 
model, we expect that data constraints on fi e in the other 
classes will be weaker than in this model (such as from 
[22] using current CMB data) , allowing for nonnegligible 
persistence of dark energy (see [4] for further demonstra- 
tions of this). Our general approach, however, does not 
rely on assuming the form for the new component or ex- 
pansion deviation. 

Figure 3 is in a sense the textbook picture of what 
Planck CMB data will say in a model independent man- 
ner about early universe expansion. For postrecombina- 
tion epochs this will improve with further ground based 
polarization measurements (especially of CMB lensing) 
and inclusion of growth of structure data. Understanding 
early expansion is in fact crucial for accurate interpreta- 
tion of large scale structure, and feeds directly into the 
early time gravitational growth calibration parameter g* 
[23]; ignorance of this can bias cosmological parameter 
estimation and tests of gravity. 

Expansion history is not the whole story as the effec- 
tive fluid behind the expansion deviations has pertur- 
bations and can have internal degrees of freedom. We 
treat the perturbations consistently - the dark radiation 
and barotropic aether models for example have sound 
speeds different from the speed of light. We do not in- 
clude viscosity, however, as the data has poor leverage on 
this [4, 24, 25]. Another difficulty for model independent 
analysis is having S < 0, since perturbations are diffi- 
cult to treat when the effective density deviation passes 
through zero; models such as nonthermal neutrinos, with 
energy densities below the standard, could realize such a 
condition. We will consider such cases in future work. 

Principal component analysis provides a valuable guide 
to the key epochs of sensitivity and the amount of infor- 
mation contributed from different times. However, we 
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FIG. 6. Three models of early dark energy are analyzed - dark radiation (top row), aether (middle row), and Doran- Robbers 
(bottom row) - with different recombination era behaviors. The leftmost column shows 8(a) built up out of PCs, with the thick 
red line giving the exact model. The second column gives the PC amplitudes m; and the third column shows the cumulative 
S/N, summing Eq. (17) over the first i PCs in quadrature. The rightmost column shows the bias squared (falling curve), 
variance (rising curve), and risk squared (top curve) when including the first i PCs. One might choose to keep those PCs that 
either give the largest jumps in S/N, or all those up to the minimum in the risk curve. 



emphasize and demonstrate that the raw uncertainty <jj 
on an eigenmode has very limited meaning; the first 15 
modes ordered by <7j can give a highly inaccurate recon- 
struction relative to a smaller number of modes ordered 
by signal to noise. Redshift bins can be more clearly in- 
terpreted. Employing the best aspects of each can result 
in physically clear, well characterized expansion history 
constraints. 
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FIG. 8. PCA of three models of early dark energy - dark radiation (top row), aether (middle row), and Doran- Robbers (bottom 
row) - is compared to the exact models in terms of the observable CMB power spectrum. The left column shows the deviation 
in the temperature power spectrum for the sum of the first 4, 8, 12, 15 modes using S/N ordering of the modes (see Fig. 6). 
The x 2 °f the deviations summing over multipoles is in the middle panels, with the dashed curve using ordering by Oi instead. 
Reconstruction of the theoretical 8(a) appears in the right panels; note how only certain epochs need be fit well to reproduce 
the observable CMB. For the aether model (middle row) dashed curves show also the cr^ ordering results for the N = 8, 15 
cases, which do much worse than S/N ordering. 



